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Abstract 

We consider nonlinear viscoelastic materials of differential type and 
for some special models we derive exact solutions of initial boundary 
value problems. These exact solutions are used to investigate the 
reasons of non-existence of global solutions for such equations. 
Keywords: Viscoelasticity of differential type; nonlinear viscosity func- 
tions, exact solutions, existence and uniqueness of boundary value 
problems. 

1 Introduction and basic equations 

In empirical rheology the standard sketch of the one dimensional KelvinVoigt 
model of linear viscoelasticity is of a purely viscous damper and a purely 
elastic spring connected in parallel. In continuum mechanics, this model may 
be recast in the class of the simple materials denoted as viscoelastic solids of 
differential type. For such materials, in the isotropic case, the Cauchy stress 
tensor, T, is defined by the constitutive equation 

T = T{B;D), (1.1) 

where B = FF^ and D = {FF^^ + F-^F^)/2 are the left Cauchy- 
Green strain tensor and the stretching tensor respectively (Truesdell and 
Noll, 1965). (Here F is the gradient of deformation). Therefore B stands 
for the elastic spring and D stands for the viscous damper. 

Nowadays there is an increasing interest in the study of nonlinear models 
of viscoelasticity and this because in the experimental studies of rubber-like 
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materials and soft tissues the viscous and nonlinear effects are prominent 
and clearly measurable (Lakes, 2009). For these reasons the possibility to 
elaborate feasible and reliable mathematical models of nonlinear viscoelas- 
ticity is a major problem of contemporary continuum mechanics (Wineman, 
2009). In this framework the class of constitutive equations (11. ip seems to 
have a central role in view of its direct and clear connection with the linear 
Kelvin- Voigt model. 

The aim of the present note is to show that many of the nonlinear models 
contained in the constitutive equation (II. ip are not acceptable because they 
are clearly unphysical. 

To make the algebra simple as possible we shall restrict our attention to 
the following family of shearing motions 

x = X + u{Z,t), y = Y, z = Z. (1.2) 

Moreover, it is convenient to consider only incompressible materials. In the 
case of compressible materials some complications arise, because the family of 
motions (11.21) is in general not admissible (Destrade and Saccomandi, 2005) 
and because the constitutive issues are more complex. 

In this framework, it is straightforward to check that the balance equa- 
tions (in the absence of body forces) reduce to the scalar equations 

putt = dTu/dZ + Pa;, Py = 0, Pz = 0, (1.3) 

where p is the (constant) density of the material and p is the Lagrange 
multiplier associated with the constraint of incompressibility {detF = 1). 
For an isotropic material the classical representation formula (Truesdell and 
Noll, 1965) gives 

Ti3 = g(A'') + iy{K\ Kl KKt)Kt (1.4) 

where K{Z,t) = uz is the amount of shear, Q = Q{K'^) is the generalized 
shear modulus and u = h'{K'^, Kf, KKt) is the viscosity function. The un- 
knowns of our problem are the pressure field p = p{x, y, z, t) and the shearing 
motion u = u{Z,t) (or equivalently the amount of shear K). 
From the three scalar equations (II. 3p we first deduce that 

p{x,t) = f{t)x + g{t), 

where g{t) is an unessential function. Then by derivation, with respect of Z, 
of the first equation in (11.30 we obtain 

P^" = 7^ [Q{K^)K + y{K\ Kl KKt)Kt] . (1.5) 
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Once that fll.51) has been solved for K{Z, t) and we integrate to obtain t) 
from (II .31) it is possible to recover the pressure field, i.e. the function f{t). 

Assuming that there is a characteristic time, w"^, such that r = cut, a 
characteristic length, L ,such that C, = Z/L and assuming that the limit 

/io = lirriK^oQ 

is a positive number, we define the quantities 

6 = , Q = —, y= — , 

such that the dimensionless version of (11. 5p reads 

d 



Q{K')K + i){K\io'K^,ioKKr)Kr ■ (1-6) 



When Q = 1, u = uq choosing u = po/i'o, (i.e. i> = 1) from (II. 6p we 
recover the classical linear Kelvin- Voigt material 

eKrr = K^C + Kr((.. (1.7) 

If in (ll.7p we let e — ?■ we obtain the quasistatic approximation. This 
approximation is fundamental to model the usual creep and recovery exper- 
imental tests (Lakes, 2009). In Pucci and Saccomandi (2010) it has been 
shown that when: 

• the empirical inequalities are in forc^ (Beatty, 1987); 

• the infinitesimal shear modulus /iq is positive and finite; 

• the viscosity function is such that 

lim u = uq, 

where z/q is a positive constant; 

then the constitutive equations (II. ip are able to describe, at least qualita- 
tively, the mechanical behavior observed in standard quasi-static experimen- 
tal tests. 



""^In this case the generahzed shear modulus, Q, is positive for any value of the amount 
of shear and the strain-energy density function associated with the elastic part of equation 
(jl.6p must be convex. 
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The quasi-static approximation is a singular perturbation of the full dy- 
namical equation fll.6p and therefore is an approximation valid for large times. 
For this reason the global existence of solutions of f ll.6p is a main issue. The 
existence, uniqueness and regularity of solutions for initial-boundary value 
problems (IB VPs) associated to (11. 5p (or its dimensionless version (II. 6p ) have 
been studied by several authors. A state of art of the various results both 
for classical and weak solutions is provided in (Tvedt, 2008). 

The main problem with (II. 5p is that the existence of smooth solutions 
global in time is not guaranteed when the viscosity function is nonlinear. 
Indeed Tvedt points out that: 

To my knowledge, all existence results for weak solutions of any 
kind to (1.1) in multiple space dimensions, including the Young- 
measure valued constructed in [6], require at least 

\S(p,q)\<C{l + \p\ + \q\) 

for some constant C. This condition is difficult to justify on phys- 
ical grounds. 

(Using the present notation 5* = T13, p = K, q = Kt.). 

It seems that the consequences of this negative result quoted by Tvedt 
are not full appreciated by workers in nonlinear viscoelasticity. For example, 
if we have not a global existence result it is clear that we can consider a 
quasi-static approximation. Our show note is centered on the impact of this 
result in very simple problems of viscoelasticity. 

The plan of the paper is the following. In the next Section we derive 
some simple exact solutions for simple IBVPs with very smooth initial and 
boundary conditions of special cases of (II. 6p . For all these solutions there is 
a finite time t* associated with the blow-up or the extinction of the solution. 
In Section 3, we suggest that the catastrophic problem associated with t* 
is located in the initial layer of the solution (i.e. for short times). By con- 
sidering the case where the viscosity function depends only on the amount 
of shear we prove that the dynamics in the initial layer is governed by an 
equation equivalent to a nonlinear diffusion equation. The blow-up charac- 
teristics of such an equation have been studied into details in (Galaktionov 
and Vazquez, 1995) and the conditions derived there show the non-existence 
of global solutions. The last Section is devoted to concluding remarks. 
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2 Some exact solutions 

We start by considering the problem associated with some viscous function 
that does not have an upper bound. To this end we consider the following 
three nonlinear polynomial viscous functions in turn 

1),{K^) = 1 + Z>eli^^ v,{Kl) = 1 + v.^Kl, v^{KKr) = 1 + VmlKKr. 

For all these viscosity functions we have that \im.x^Q,Kr^o P = 1 > but 
none of these functions are bounded from above. 
We start considering, for G [0, 1], the equation 

eJC. = /^cc + Krcc + ^ei [K^K,]^^ , (2.1) 

i.e. an equation corresponding to a material whose elastic part is neo- 
Hookean and with a super linear viscosity function in the amount of shear 
defined as i>e- 

This equation has an interesting and simple exact solution that may be 
computed as 

jaC,r) = a(r)C + /3(r), (2.2) 

where 

da ^ s df3 ^ On 

e— = ^Ueia + ci, e— = 2z/eia p + C2, 
or uT 

(here Ci and C2 are two constants of integration). If we consider Ci = C2 = 

and 

2 2 Y e — I'eiT 

we obtain 



^(C, r) = Ll^—iC -b^ HC r) = \J^—ie - C). (2.3) 
2 y e — z/giT 2 4 y e — z/gir 

(Here m = u/L). 

The shearing motion in 02.31) satisfies the following set of smooth initial 
and boundary conditions 

n(0, r) = n(l, r) = 0, n(C, 0) = (C' - C)/4, u^C, 0) = i>i(C' - C)/(8e). 

In this way we have provided a simple example of non-existence (because of 
finite time blow-up) of the solution for a very regular IVBPs for equation 
dSD in [0,1]. 

If now we consider the equation 

eKrr = i^cc + + hi (K^) , (2.4) 
5 



a similar solution may be obtained using the ansatz in (I2.2p for the amount 
of shear. The exact solution is now 

This solution may be used to solve the boundary-initial value problem 

u{0, r) = u{l, r) = 0, u{C, 0) = C - C, «.(C, 0) = 'T^iC - 0- 

Here we record that the reason why the existence of the solution fails in finite 
time is the extinction phenomena. 
The last equation we consider is 

eK^r = K^^ + ir^cc + {KKl)^^ . (2.6) 

Here the mathematics is more involved. Let us consider that the shearing 
motion is in the form 

n(C,r) = a(r)(C^-C). 

Then the time dependent function a is determined by the ordinary differential 
equation 

-^a" = 2Aaa'\ 

The solution of this second order ordinary differential equation is given by 

Erf (a) = ciT + C2, 

where 

a = 2 Y «, 

Erf is the standard error function and Ci, C2 two integration constants. This 
exact solution solves the following IBVP (in the normalized variables) 

«(0,r) = n(l,r) = 0, S(C, 0) = 5(C' - C), 

with < 5 < 1 and £t^(C,0) = /i(C^ — C) "where /i is an ad hoc constant 
determined by a' . From 

Erf(a) -Erf(5) = ciT, (2.7) 

it is clear that this solution blows up in a finite time. 

To summarize, the following remarks apply to all the three examples we 
have proposed. First of all in the case e — )■ all the exact solutions we propose 



6 



degenerate and are of no interest. Second, for all three solutions there exists 
a critical time r* such that for r > r* the solution does not exists anymore. 
For the solution in (12.31) and (12.71) we record a blow-up phenomenon and for 
the solution (12. 5p we have an extinction phenomenon. 

Therefore we have provided explicit examples of the lack of global exis- 
tence for equations with nonlinear viscous function and a neo-Hookean elastic 
term. For all these solutions the inertia term (i.e. e 7^ 0) is crucial in delin- 
eating why the global existence fails. 



3 The initial layer 

We point out once again that in the quasi-static limit all the models we have 
considered in the previous Section are "good" models. This fact suggests 
that because the quasi-static limit is the result of a singular perturbation 
procedure, the failure of well-posedness must occur for short times or, using 
the terminology of singular perturbations in the initial layer. 

When we consider a general viscosity function defined as P = v{K'^) and 
we introduce the stretched time f = r/e, it is possible to obtain easily, for 
any regular viscosity function, the equation governing the dynamics in the 
initial layer. Restricting the attention to this class of materials it is then 
possible to rewrite (II. 6p as 

Kff = ^ [eQiK')K + i){K')K,j , (3.1) 

and when e — )■ we derive the equation 

where U = J u{K'^)dK. 

By a simple integration the above equation is rewritten as 

Kf = ^ + ^{C). (3.3) 

Equation (13. 3p is a nonlinear diffusion equation. This equation has been 
considered by many authors and the necessary and sufficient conditions for 
the complete blow-up and extinction of this equation have been considered in 
(Galaktionov and Vazquez, 1995). It is interesting to note that for equation 
(13. 3p we have blow-up of the solution in finite time if the integral 
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does not converge. (Here a dash denotes the derivative with respect the 
amount of shear). This is exactly the case of the examples in the previous 
Section. 

It is clear that there are several choices of the viscosity functions for which 
the above integral may converge, but the usual higher order polynomial mod- 
els that are used in the standard fitting procedures do blow up. We remark 
that a model with a constant viscosity, in the initial layer, is associated with 
a linear heat equation this explain why it is possible to find a unique global 
solution for such a model in the class (11. ip . 

We remark that via a multiple scale methods is possible to deduce the 
equation in the initial layer for more general viscosity functions, but the case 
we have just examined is sufficient to point out the problems inherent with 
the constitutive equation (II. ip when the sub linear growth condition pointed 
out by Tvedt (2009) is not satisfied. 

4 Concluding remarks 

If we believe that materials of differential type have a "good" quasi-static 
behavior, the fact that a global solution for the full dynamical equations 
does not exist is a major problem. Indeed, if a solution blows up, or is 
associated with an extinction phenomena, in finite time it is not possible to 
consider the asymptotic limit for t — >■ oo and therefore it is not possible to 
consider the quasi-static approximation. 

Clearly such situation may be arranged considering a sort of "cut-ofF' 
for the usual polynomial constitutive equations. For example we may be 
consider a constitutive equation such as 



We cannot be satisfied with this approach for two reasons. First of all we 
introduce some oddities in the quasi-static regime as it is possible to check by 
analyzing the corresponding creep problem that loose regularity. Moreover, 
we do not solve the "unphysical" situation generated by the fact: greater 
viscous effects lead to greater problems with the regularity of the solution. 

This situation seems to suggest that there is something missing in the 
"dynamics" of the nonlinear viscoelastic models we are considering. It seems 
that the model (II. ip must be refined taking into account phenomena that 
at first sight seemed to be irrelevant but indeed are fundamental. These 
phenomena seems to be of the same order of the inertia terms and therefore 
disappear when e — )■ 0. 




1 + u^iK^ if < K*; 
V, if K"^ > K*. 
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A possibility is to consider the presence of objective derivatives of the 
stretching tensor D i.e. higher order models of differential type where the 
constitutive equations not only depends on the stretching tensor D but for 
example also on the subsequent Rivlin-Ericksen tensors (Truesdell and Noll, 
1965). These tensors may introduce some additional terms that regularize 
the equations and allows to avoid the problems presented here. This direc- 
tion is outlined in the context of non dissipative materials in (Destrade and 
Saccomandi, 2006). 
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